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Abstract 



In this paper we develop a polymer expansion with large/small field conditions for the mean 
J5 ' resolvent of a weakly disordered system. Then we show that we can apply our result to a two- 

dimensional model, for energies outside the unperturbed spectrum or in the free spectrum provided 
the potential has an infra-red cut-off. This leads to an asymptotic expansion for the density of 
^ , states. We believe this is an important first step towards a rigorous analysis of the density of states 

Q ' in the free spectrum of a random Schrodinger operator at weak disorder. 



1 Introduction 



^~~^ ] In the one-body approximation, the study of disordered systems amounts to the study of random 

I ■ Schrodinger operators of the form 
(N H = Ho + XV (1) 

o ■ 

t~^ ' where Hq is a kinetic term (i.e. a self-adjoint or essentially self-adjoint operator corresponding to 
^» ■ some dispersion relation, typically a regularized version of —A) and F is a real random potential (in 
"^ , the simplest case, F is a white noise). We work on a ultra-violet regular subspace of £^(IR ) and we 
ri I restrict ourselves to A small so as to see XV as a kind of perturbation of the free Hamiltonian. 

The properties of H are usually established through the behavior of the kernel of the resolvent 
operator or Green's function ( [ pi'hol , [ Aiz |, Era ]) 



o 

p.: GE{x,y)=<x\j^^\y> (2) 



k>^ , For instancej the density of states is given by 



p{E) = - lim Im GE+ieix, x) (3) 

vr e^O 

The important point is that, in the thermodynamic limit, the system is self- averaging, i.e. mean 
properties are often almost sure ones. Thus the problem can be seen as a statistical field theory with 
respect to the random field V. In Statistical Mechanics, functional integrals in the weakly coupled 
regime are controled through a cluster expansion (or polymer expansion) with small field versus large 
field conditions, the problem being then to control a Boltzmann weight ( [ [Bry| ] , | Rivl| ). 



In the first part of this paper, we derive a resolvent cluster expansion with large field versus small 
field conditions assuming that V satisfies some large deviation estimates. This would allow to prove 



the existence and the regularity of the mean Green's function (theorem |l|) and to get an asymptotic 
expansion for the density of states. 

In the second part, we show that the hypothesis of theorem |l] are satisfied in the case of a 2 
dimensional model with a rotation invariant dispersion relation and an infra-red cut-off on the po- 
tential. From the point of view of Renormalization Group analysis, our results allows to control the 
model away from the singularity, i.e. to perform the first renormalization group steps and therefore 
to generate a fraction of the expected "mass" . 

2 Model and results 

2.1 The model 

In IR we consider 

H = Ho + XV (4) 

where ^ is a gaussian random field with covariance ^ whose smooth translation invariant kernel is 
rapidly decaying (we will note the associated measure d/x^). Because ^ is smooth, dfi^ as a measure on 
tempered distributions is in fact supported on C°° functions. We suppose also that Hq has compact 
support so that we do not have to deal with ultra-violet problems. We construct the finite volume 
model in IR / . -_rf by replacing ^ and Hq by their "A-periodization" 



Ad 



pe^2'* zeAZ"* 



H^o\x,v) = ...= Y. Ho{x-y + z) (6) 



Then we define 

GA,eiE,X,V) = -^ — (7) 

H^o> + XV-{E + ie) 

GA,e{E,X) = Jdfl^^{V)GA,e{E,X,V) (8) 

where d//^^ can be considered either as a measure on C°° ( IR /. -^d 1 or as a measure on C°°(IR ) which 
is supported by the space of A-periodic functions. In the same way, GA,e will be considered as an 
operator either on >C^(IR / ^-^d) or on i2^(IR ). One can note that in momentum space, because of 
the cut-off, the problem reduces to a finite dimensional one. 

Because V is almost surely regular, its operator norm as a multiplicative operator is equal to its 
C°° norm (it is easy to see that ||F|| ^ ||^||ooi equality can be obtained by taking test functions fn,x 
such that fnx^ ^x)- Therefore V is bounded and self-adjoint. Then GA,e{E,X,V) is almost surely 
an analytic operator-valued function of A in a small domain (depending on V) around the origin. This 
domain can be extended to a ^-dependent neighborhood of the real axis thanks to the identity (for 
|A — /x| small enough) 

Ga,,{e, /x, V) = GaAe, x,v)h+Y{x- fiT [vGaAe, a, V)A (9) 

In the same way, GA,eiE, A, V) is analytic in E. One can also check that GA,e{E, A, V) has a smooth 
kernel and is integrable with respect to d^^^^. Furthermore, GA,e{E, A) will have a translation invariant 
kernel because d^^^ ^^ translation invariant. 



2.2 Main result 

We introduce a function 9 which satisfies 

• ^ is an odd C°° function, increasing and bounded 

• for any x, |^(x)| ^ \x\ 

• for any |x| < 1, 6{x) = x 

• the £°° norm of its derivatives does not grow too fast 

Then for ;U > 0, we define the operators Ca,/^, ^A,/i and U\^^ through the Fourier transform of 
their kernel 

CkI{p) = Hi^\p)-E-^^i (10) 

^^Ap) = — — ^ -2 (11) 

U]^M = Dl,iP)CA>) (12) 

Given any characteristic length L we can divide the space into cubes A of side L and construct an 
associated C^ partition of unity 

1 = EXA (13) 

A 

where XA has support in a close neighborhood of the cube A {e.g. on A and its nearest neighbors). 
This decomposition induces an orthogonal decomposition of V into a sum of fields Va with covariance 

C^{x,y) = Jdze/\x - z)xA{z)e/\z - y) (14) 

For simplicity we will pretend that ^ and ^^'"^ have compact support, so that Va is almost surely 
supported on a close neighborhood of A, moreover we will take that it is restricted to A and its nearest 
neighbors. The generalization to a fast decaying ^ can be easily obtained by decomposing each Va 
over the various cubes and write more complicated small/large field conditions that test the size of 
V^ in the various cubes. This leads to lengthy expressions that we want to avoid. 

Finally, we note d\ the distance in IR / . -.^ 

d\{x.,y) = min \x — y + z\ (15) 

In the following, C or 0(1) will stand as generic names for constants in order to avoid keeping 
track of the numerous constants that will appear. Furthermore we will not always make the distinction 
between a function and its Fourier transform but we will use x, y and z as space variables and p and 
q as momentum variables. 

Theorem 1 

Suppose that 



• ^ is smooth and has fast decay 

• Q = sup— -T / il^{x,y)dxdy exists 



for all E G [£^i, £'2] oind all n, C^, I?^ and U^ have smooth kernels with fast decay over a length 
scale L. 



for all 71-1, we have Cm such that for all A and all triplets (Ai, A2, A3) 

\\xa.da,,va,da,,x^,\\ ^ |rTi^H^(Ar;A^rMrTL^H:(A^;A^ ^^^^ 

there are constants Cq, Ci, k > and a > such that 



VA s; 00, Va > 1, VA, Pa (||Z?a,^FaI?a,^|| ^ aCo) ^ de"^'^ ^" (17) 

where -/Pa(-) denote the probability with respect to the measure djj,^^ = (ddiicA f^oo = C) 

IPa{X) = Jdii^^iV) lx{V) = Ma(^) (18) 

Then let fiQ = L~'^''^C/ , jjl = XfiQ and 

TA,e = D]^lfiA,sD]^^ (19) 

For all A ^ Aq = 0(1) and for all e small enough (in a \-dependent way), TA,e(-E, A) is uniformly 
bounded in A and admits the following development (in the operator norm sense) 

ln^TA,e{E, A)lf,^ = ln^T{E + ie, A)!^^ + O (^j^ (20) 

where Qa = [— A-*^'^; A^'^]"^, and 1q^ is the characteristic function of Qa- 
Furthermore we have the following properties 

• T has a smooth, translation invariant kernel 

• Ta^e o.iT'd T have high power decay 

3no large, BCrino) such that V(A, A'), ||1aTa,,1a'|| ^ ^ , r ^Z tl aai"o (21) 

[l + L--^(iA(A,A')J 

and a similar relation for T with d^ being replaced by d. 

• T[E, A) is an analytic operator valued function of E for all E in ]Ei, E2[ with a small X-dependent 
radius of analyticity. 

• T{E, A) is a C°° operator-valued function of A and admits an asymptotic expansion to all orders 
in A. which is the formal perturbative expansion of 

dfi.{V) e^<i'«"'i> e^^o<^'«"'i> -i — (22) 

(<> denotes the scalar product, i.e. < f,Af >= J f{x)A{x,y)f{y) dxdy) 

This theorem is formulated in a rather general way so as to apply with minimum transformation to 
various situations (lattice or continuous models) and in any dimension. Then we construct a concrete 



example with a two-dimensional model. One can also refer to |MPR| for a d = 3 case. 



2.3 Anderson model with an infra-red cut-off in dimension d=2 

We consider 

H = -A^ + XmVm (23) 

where 

• Az"^ is a ultra-violet regularized inverse Laplacian, i.e. there is a Cq° function rjjjy equal to 1 on 
"low" momenta such that 

A-\p) = ^ (24) 

We will note p^ instead of — A^, the UV-cutoff being then implicit. 

• we are interested in the mean Green's function for an energy E = 0(1) 

• r]E is an infra-red cut-off which enforces 

|p2_E| ^ AA^jlogAp (25) 

for some large constant A 

• V has covariance E, which is a Cq° approximation of a (5-function 

This corresponds to the model away from the singularity p^ = E in a multi-scale renormalization 
group analysis, we will show that it generates a small fraction of the expected imaginary part which 
is 0(A2). 

Let M^'^ be an even integer greater than 2, we define jo G IN such that 

M-^o :^ inf |p2 - E\^ M-(J"-i) (26) 

Supp(??b) 

Next, we construct a smooth partition of unity into cubes of side M^° (they form a lattice SDj^) 
and we construct the fields Va's accordingly. 

Theorem 2 

There exist constants Cq and Ci such that for any A, a ^ 1 and A G JDj^ we have 

Pa {\\DA,^.VEVAVEDAJ\ ^ aCoJo^P'^/^) ^ Cie-^y^'''^' (27) 

Furthermore theorem Q| applies and Ge is asymptotic to its perturbative expansion so that it behaves 
more or less like 

^^ ^ p^-E-ir]EO{X^\\ogX\-'^)r]E ^^^^ 

It is easy to extend this result to the case of a rotation invariant dispersion relation and for energies 
outside the free spectrum not too close to the band edge. In this case, the cut-off is no longer needed 
so that the result apply to the full model. 



3 Resolvent polymer expansion with large field versus small field 
conditions 

3.1 Sketch of proof for theorem |I] 

We give here the global strategy for proving theorem |l|, the main ingredient being the polymer expan- 
sion that we will detail in the following. 

First we recall (without proving them) some quite standard properties of gaussian measures. 

Lemma 1 Complex translation 

Let X be a gaussian random field with covariance C and let dfic be the associated measure. For any 

regular functional J-{X) and any function f G Ran C , we have the following identity 

JdficiX) HX) = e5</'^"^/> jd^ic{X) T{X - if) e*<^.^- V> (39) 

Lemma 2 Integration by part 

With the same notations than above we have 

JdficiX) X{x)T{X) = jdy C{x, y) J dfic{X)jj^T{X) (30) 

Those lemmas could for instance be easily proved for polynomial functionals and extended through 
a density argument to a wide class of functionals. ■ 

Our starting point is obtained by applying lemma || with / = /iol. 

G^,,iE + z,X) = /d/ie,(y)e#<i'«A'i>+^'^o<^'«A^^>— ^^ ^ (31) 

J H^^> -{E + ifi) + XV -ie-z 



T^AE + z,X) = /d^g,(y)e^<^-^A^^>+'^o<^-^A^^> ' ---—^(32) 



On one hand we earned something because now the resolvent operator in the integral is bounded 
in norm independently of e (in the following we will note z instead oi z + ie and show convergence for 
any z such that \z\ <C ^, this would allow to prove analyticity in z). But on the other hand we have 
a huge normalization factor to pay. However, we can remark that this normalization factor is in fact 
equivalent to a factor e per L-cube. 

Most of the demonstration amounts to a polymer expansion of Ta^e, i.e. we write Ta,£ as a sum 
over polymers of polymer activities 

Tout,in = XAout'^A,eXAi„ (33) 



-'-out,in — XAou 



[l + L ^(iA(Aj„,Ao«i)J"-" ,j^ 



XA,„ (34) 



where ci and C2 are small constants, Ty has decay in the spatial extension of Y and || TJ 7'(^)|| is 

Y&A 
bounded. Furthermore, G{Y) is given by a functional integration over fields Va's corresponding to 

cubes in the support of the polymer Y . This show that Ta,£ is bounded and has a high power decay 
uniformly in A. 

Next, when we consider InA^A.elfiA ^^ ^^^ divide the sum over polymers into a sum over polymers 
with a large spatial extension (say A^'^) and sum over "small" polymers. The large polymers will 



have a total contribution small as A^^ to some large power. For the small polymers, since we are far 
away from the boundaries, their contribution calculated with d^i^^ will be equal to their contribution 
calculated with d/i^ up to a factor A~". In this way we can prove the development (|20|). Smoothness 
of the kernel will be obtained because we will show that we can write 

T{Y) = UA,^f{Y)UA,,, (35) 

The convergence for any z <^ fi allows to show analyticity (we write z-derivatives as Cauchy integrals 
so that we can show that they all exists and do not grow too fast). Then an asymptotic expansion 
can be generated through the repeated use of resolvent identity. 
Finally, for the density of states, we just need to remark that 

G(0, 0) = [dp dq G{p, q)=<6,G6> (36) 

where (^ is a regularized 5-function because of the presence of the ultra-violet cut-off. Thus an asymp- 
totic expansion for G with respect to the operator norm will yield an asymptotic expansion for the 
density of states. 

3.2 Improved polymer expansions 

Cluster expansions in constructive field theory lay heavily on a clever application of the Taylor formula 
with integral remainder. Writing the full Taylor series would amount to completely expand the 
perturbation series, which most often diverges, and therefore should be avoided. A rather instructive 
example of minimal convergent expansion is the Brydges-Kennedy forest formula: you have a function 
defined on a set of links between pair of cubes and you expand it not on all possible graphs but only 



on forests (cf | Bry |). 



For more complex objects a way to generalize such a formula can be found in |AR], and we refer 
the reader to it for a more careful treatment and for various proofs. Let us assume that we have a 
set of objects that we call monomers. A sequence of monomers will be called a polymer, then we will 
expand a function defined on a set of monomers into a sum over allowed polymers. 

To be more precise, let ^ be a set of monomers, we define the set y of polymers on X as the 
set of all finite sequences (possibly empty) of elements of X. Then a monomer can be identified to a 
polymer of length 1. The empty sequence or empty polymer will be noted 0. We define on 3^ 

• a concatenation operator: for Y = {Xi, . . . , Xn) and Y' = {X[, . . . , X'^,), we define 

YUY' = {X,,...,Xn,X[,...,X:,,) (37) 

• the notion of starting sequence: we say that Yi is a starting sequence of Y (equivalently that Y 
is a continuation of Yi ) and we note Yi C Y iff there exists Y2 such that Y = Yi U I2 

Then we call allowed set (of polymers) any finite subset A C y such that 

• VY, Y' Y' C Y and Y G ^ ^ Y' G ^ 

• yX,Y,Y' YcY'andYUX^A=^Y'uX^A 

the first condition implies that 9 € A whenever A is non-empty. Finally, for Y belonging to some 
allowed set A, a monomer X is said to be admissible for Y (according to ^) iff Y U X G ^. 



Lemma 3 

Let X = {X} he a set of N monomers and y the set of polymers on X. We assume that we have an 
indexation of IR by X, i.e. a bijection from X to {1, . . . , A^} so that an element of IR can be noted 

z = {zx)x£X- 

For T a regular function from IR to some Banach space B and an allowed set Ady , the polymer 
expansion of T according to A is given through the following identity 



^(1) ^^(1,...,1) 



E E 

n^O Y=iXi,...,Xn)eA 



dhi . . . dh„ 

l>hi>...h„>0 



where z{Y, {hi}) is given by 



n 

\X£Y 



d 

dzx 



:F[z{Y,{h,})] 



(38) 



zxiY,{hi}) 



if X is admissible for Y 

1 if X is not admissible for 

hi if X not admissible for Y and X = Xj for some j, 

in which case i = m.ax{j/X = Xj} 
hi with i = nim{j/X not admissible for {Xi, . . . Xj)}, otherwise 



(39) 



Proof 

The proof is made through an inductive iteration of a first order Taylor formula. We start with 
J-"(l) and put a common interpolating parameter hi on all admissible monomers for the empty set, 
i.e. we make a first order Taylor expansion with integral remainder of .F hizi + {1 — zi) between 
and 1, with zi being the vector with entries 1 or according to whether the corresponding monomer 
is admissible or not. Then each partial derivative acting on J^ can be seen as taking down the 
corresponding monomer so that terms can be seen as growing polymers. The iteration goes as follow: 
for a term of order n corresponding to a given polymer Y and having n interpolating parameters 
1 > /iQ > • • • > ^n > we put a common parameter /i„+i interpolating between and /i„ on all 
monomers admissible for Y. It is easy to check that the process is finite since A is finite and that one 
obtains the desired formula. ■ 

In the following our monomers are sets of cubes (that we call the support of the monomer) and 
links between those cubes. When we take down a polymer, we connect all the cubes in its support 
and maybe some more cubes. Thus a polymer is made of several connected regions, we will say that 
it is connected if it has a single connected component. The rules of admissibility will be to never take 
down a monomer whose support is totally contained in a connected region. 

In this case, one can show that the interpolating parameters depend only of the connected compo- 
nent to which the corresponding monomer belongs so that one can think to "factorize" the connected 
components. We define generalized polymers as sets of connected polymers. Then a generalized poly- 
mer Y = {Yi, . . . , y^} is allowed if the polymer Yi U . . . U 1^ is allowed (this does not depend of the 
order of the Y^s). Equation (pi) becomes 



^(1) 




dhj . . . dh'l^ 

l>h^>...h"^>0 




\xeY 



^{Y,{hl}) 



(40) 



where the sum extends on all allowed generalized polymers, and z{Y,{h-'-}) is given by 



^x(Y,{H}) 



' if X is admissible for Y, i.e. for Yi U . . . U 1^ 
1 if X is not admissible for 

hj li X = Xj for some i and j (41) 

hi where X is not admissible for Yi and 

j = minjfc/X not admissible for {Xj, . . . ,X^)}, otherwise 



3.3 Large/small field decomposition 

Semi-perturbative expansion (like cluster expansions) are convergent only when the "perturbation" is 
small (in our case the operators Va's). Thus it is very important to distinguish between the so called 
small field regions where perturbations will work and the large field regions where we must find other 
estimates (they will come mostly from the exponentially small probabilistic factor attached to those 
regions). 

We take a Cq^ function e such that 

• Os^e^ 1 

• Supp(£) C [0, 2] 

• "l[o,il = 1 



Then for each A we define 



l-^A.M^A^^.MI 



eAiVA)=e r ;^X-i/4^^ I and 7?a = 1 - ^A (42) 

where a = 0(1)- Then we can expand 

A N^O n={Ai,...,A]v} \Aen I \A^n / 

where Vl is the large field region whose contribution will be isolated through the following lemma. 

Lemma 4 

Let Q be a large field region made of N cubes Ai, . . . , A^v and A any operator such that 

yn c{l,...N}, ^ + E ^i «s invertible (44) 

ieD 

(Bi stands for B^- = ADa^^Va^-^a.^J. 
We have the following identity 



H n=0 iie{l...N} i2e{i...Af} i„G{i...iV} 

i2^{n} »n^{n-..in-l} 



where 



0, = B,-\ E bA \^ B, (46) 

ViG{l...p} ) ^^ 1^ ^^ 

ie{i...p} 



Proof 

The proof relies on resolvent expansion identities 



1 



"l-B- ' 



I 



A + B A V A + B 

We show by induction that for all m G {1, . . . , A^} we have 

1 m— 1 

A "A 



' y 



A + B y A 



A + Y.Bi 

Ihrn. 



n=0 (ii,...i„) 



A 



E 



1 



(n,---im) 



A+Ei?;^*"^^™"'^"""^'A 



The case m = 1 is obtained by a resolvent expansion 



1 



1 



E 



1 



-5, 



1 



^ + ESi A ^ A + ZBi ''A 
Then we go from m to m + 1 with 2 steps of resolvent expansion. We write 

1 f . ^r^ 1 „ \ 1 



A + EBi 



I 



E 



-Bi 



A -I- \^ R- ^"^+1 

fe = l 



/ 



E 



1 



«m+l^{n---Jm} 



A + E^i 



-B 



«m+l 



.4 



^-E^^.- 



v 



=^ ^ + E^. 



Z=l 



Finally, for m = N we make a last resolvent expansion on the rest term Ri\f by writing 



1 



1 



A + J:B^ A 



tU-E^^ 



1 



A + EBi 



(47) 

(48) 
(49) 



(50) 



(51) 



(52) 



(53) 



If we look at 

XAout ^^y^ XAi„ (54) 

and fix {Aj^ , . . . Aj„}, we can see that summing over the sequences (ii, . . . i„) and choosing a particular 
term for each Op amounts to construct a tree on {Ai„, Aout, \i, ■ ■ ■ Aj„}. 

We define an oriented link I as a couple of cubes that we note (l.yj.x), then C is the set of 
oriented links. Given two cubes Aj„ and Aout and a set of cubes 0, = {Ai, . . . A„} we construct the 
set Tr(Aj„, Aout, ^) of oriented trees going form A^ to Aout through Q as the sequences (^i, . . . In+i) G 
£"+^ which satisfy 

• h.x = Ain 

• Ui+i-y = Aout 

• VA; G {1, . . . n}, /j^.y G il 



10 



• VA: G {2, . . . n + 1}, k-x e {h.y, . . . h-i-y} 

• VA; G {2, . . . n}, h-y {h-y, ■ ■ ■ k-i-y} 

Then we have the following equivalent formulation of lemma |^. 

Lemma 5 

Let il. be a large field region made of N cubes Ai, . . . , A^v and A any operator such that 

yDc{l,...N}, A + J2 Bi is invertible (55) 

We have the following identity 

^i2_^-^i ra=0 n'cii rGT^(A,A',n') 

n'={A[,...A'„} r=(ii,...i„+i) 

XA'^On{ln+l-X,ln-y)-7 ■ ■ ■ Oi{l2.x,h-y)^XA (56) 

Op(A„ A,) = BaM.a, - Ba- ^^— i?A. (57) 

A+ y ^ Bj 

ie{i...p} 

The proof being just a rewriting of lemma ^ is quite immediate. ■ 

Thanks to this lemma we can factorize out the contribution of the large field region, then we need 
to extract spatial decay for the resolvent in the small field region. However a kind of Combes-Thomas 
estimate ( [|CT| ) would not be enough because of the normalization factor that we must pay. For 
this reason, we will make a polymer expansion to determine which region really contributes to the 
resolvent. 

3.4 Polymer expansion for the resolvent in the small field region 

For some large field region $7, we want to prove the decay of 

^ = RUa,^, (58) 



where 



^A,i + A E ^A,M^A^A,/^ - ^^A,M 



and get something to pay for the normalization factor. 

We define the set C of links as the set of pair of cubes, and £(i7) as the set of links which do not 
connect two cubes of il. Then for / = {A, A'} we define 

Ql = A(xAf^A,^i^A,At^A'^A,M + XA't^A,/i-CA,^VAl?A,/x) -2(XaCa,^XA' +XA'C'a,mXa) (59) 

Ri = xaRxa' + XA'RxA (60) 

Ui = XAUA^f^XA' + XA'Ua,^xa (61) 

with the convention that Va = if A G fi. 

For any fixed Iq = (Aq, Aq) we expand i?;^ on £(^2) with the rule that for any growing polymer 

• if we have two adjacent connected components Yi and Y2 (such that dA{Yi,Y2) = 0) we connect 
the two components 

11 



• we connect Aq (resp. Aq) to any adjacent polymer component 

This allows to take into account that the operators localized on a pair of cubes have their support 
extending to the neighboring cubes. 

Let us notice that if A and B have disjoint support, we have 

111 

(62) 



I+A+B I+AI+B 

Then it is easy to see that the expansion of Ri^ involves only totally connected polymers which 
connect Aq to Aq, because the other terms necessarily contain a product of two operators with disjoint 
supports which gives zero. We note A{Vt,, Iq) the corresponding set of polymers which is a decreasing 
function of O, i.e. 

n' cn^ A{n, lo) c A{n', /o) (63) 

Then, according to (^0|), our expansion looks like 

^ [z{Y,{h})] (64) 



^;o(^) = E E / dh,...dhJl[J-] 



/+ E ^^Q^ 

Y={xi,...,x„) xec(n) 

= E Jl[dhJll^)j-^ -^—^^--[z{Y,{h.}),0]m 

Y&Amo)^ i \x&"^xjl+ 2^ zxQx+ l^uxQx"- 

XGC(n) xeY 

Then in the second expression, we rewrite the derivatives as Cauchy integrals so that 

Rioi^) = E JlldhJuf^)-—^ ,^ ' ^v n ^''^ 

YeAinMY i Kxev^ ^^"^^x J I + 2^ zx(Y,{hi})Qx + }_^ uxQx 

xec{n) x& 

^ Y: R{Y) = I+ J2 R(Y) (67) 

YeA{nM Y£A*{n,io) 

where A*{n, Iq) = A{n, lo)/{^} 

We suppose that we fixed ni the power rate of decay in d^i^Ai, A3) of ||xAi-DA,/iVA2^A,/iXAall aiid 
IIXAiC'A.^XAall, then we have the following lemma 

Lemma 6 

For n2 = rii — 3{d + 1) and A small enough, we have 

V/o = {Ao, A'o}, Vy G ^*(0,/o) = A*ilo), 

\\y\/^ 
\\R(y)\\^T. — . 1 , ,. — ^^^r(y) mi/i y r(y)<i (68) 

" ^ ^" [l + L-idA(Ao,AoO]"2 ^ ^ yg^(^^) ^ ' ^ ' 

where \Y\ is the number of monomers in Y . 

Proof 

Since we are in the small field region 
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Then in (pq ) we can integrate each ui on a circle of radius 



(70) 



while staying in the domain of analyticity for ui and have a resolvent bounded in norm by say 2 (if A 
is small enough). Thus 



\\R{y)\\ ^ 



( 



dh\ 



dh 



€ 



l>hl>.../l|y|>0 

Al^l/4 



y\ 



n 



0(l)Ai/2 






0(l)Ai/4 



\Y\ 



[l + L-i(iA(Ao,Ao')] 



n2 



irl! 



n 



this demonstrates the first part of the lemma with 

0(l)Ai/4]'''' 
T{Y) = i 



\Y\\ 



nrx 



\X = {A^,A^} 



and 5] Tx = 0(1) 

X9A 



(71) 



(72) 



(73) 



A link / = {A, A'} can either be a true link when A 7^ A' or a tadpole when both cubes collapse. 
Our expansion rules insure that there is at most 1 tadpole per cube of SuppiY) the support of Y . If 
we forget about proximity links (that we connect also adjacent cubes) then a polymer with m true 
links and p tadpoles has a support of ?7i + 1 cubes (2 of them being Aq and Aq') and the true links 
make a tree on the support of Y . If we take into account the proximity links then two connected links 
in the tree on Y are adjacent instead of sharing a common cube, we will forget about this since it 
would induce at most a factor 0(1)' '. The links in Y are ordered, but we can take them unordered 
by eating up the 1/|^|! we have in r(y). 

Then the sum over Y can be decomposed as 

• choose m^ 1 

• choose 771 — 1 cubes {Ai, . . . , Am-i} 

• chose a tree T on {Aq, Aq, Ai, . . . , Am-i} 

• choose ^ p ^ m + 1 

• place p tad-poles on {Aq, Aq, Ai, . . . , Am-i} 

We can perform the sum on tadpole configurations because for p tadpoles, we have a factor 
0(1)A^''^ coming from the tadpoles and at most ( "^ j configurations. Thus 



H^Y) ^ E E E[i + o(i)A^/' 

Y&A*(lo) m>l{Ai,...A„_i} r 



m+1 



0(1)A 



A^ 
1/4 



\X£r / 



(74) 
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Then we fix first the form of T then we sum over the positions of Ai, . . . , Am-i- But since the cubes 
are now labeled we get (m — 1)! the desired sum. 

E W^EhDA^/^T^^T^E E (U^x) (75) 

YeA-ih) m^i ^ ^' r (Ai,...,A™._i) \xer I 

We choose Aq as the root of our tree and suppose that the position of Aq is not fixed. Then the sum 
over the position of the cubes is made starting from the leaves thanks to the decaying factors V x (cf. 



[|Riv2| ), this costs a factor 0(1)'". 

Finally, the sum over T, which is a sum over unordered trees, is performed using Cayley's theorem 
which states that there are (m + l)*"""*^ such trees. 

Y^ T{Y) < Y. [0(1)A^/^] "" ^7 +^\7' ^ ^('^^'^' ^ ' ^^^) 



YeA'Ho) m>l 



(m-1)! 



for A small enough. 

We note that we can perform the same expansion on 



1 



R' = UaI — -. ^^ ^ (77) 

UKi. + A E Da,,VaDa,, - zDl^ 

3.5 Summation and bonds on T 

We define 

Tout,in = XAout'^A,eXA,„ (78) 

We can combine equations (|3l|), (|4^), ( |5^ ) and (|67|) to write 



/0d/u^A(VA) e4<i'«A'i>+'^o<^'«A'i> J2 E ( n 



A^SsO Q={Ai,...Ajv} VAen 

E(-ir E E E E E 




\n'\=n r={ii,..J„+i) (aI',...a21) Y,6.A(n,{A«,A^}) 

f^A„„t,A^+i-RA^_^^,A-_^^^n(Wl-a;,^n-y)-RALAjif^A-,A- • • • ?1 (/2.x, ^l-y)-RA5',Af f^AJ,A,„ (79) 

where we pretend that the xa's are sharp otherwise we would have to deal with adjacent cubes but 
it's an irrelevant complication. Furthermore, for the leftmost term we made a polymer expansion of 
UA,fj.R' instead of RUh^fi so that we can write Tout, in as 



'-out An 



XAout {Ua,^ + UA,^,fUA,^?) XA,„ (80) 



The crucial point here is to notice that for any cube A, each term where A appears in fi but not 
in 0' pairs with a corresponding term where A fi and A IJSupp(l^) (i.e. A has been killed in 
every polymer expansion). Then the corresponding ea and t^a add up back to 1 so that 

Tout,in =2j(~i)" E E E E E 

n^O f7={Ai,...A„} r6rfj{A,„,A„„t,n) (Af ,...a-_^j) (A|,...A^^J {n.-.^n+i) 

r=(!i,...i„+i) (A^,...A«) " 
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VAen / \AGuSupp(yi) 

2 

The factor e"2"<^''^A i>+«mo<V',?a i> correponds to the translation of V by —i/J-o, this is equivalent 
to have translated all the Va's by —i^oXA therefore we can write it as 

2 

A 

then we can perforin the integration on all Va ^2 U (USupp(l^)) so that the normalization factor 
reduces to 

TT g^<XA,«^)-ixA>+w<VA,(Sf)-ixA> (83) 

Aef7u(uSupp{yi)) 
This amounts to pay a constant per cube of il U (IJ Supp(yi)), this is done in 0, with a fraction of 
the probabilistic factor coming from the large field condition and in IJ Supp(l^) with a fraction of the 
factor A-^ l^^l/^ coming from the R's. 

The sums over the various 1^'s are controled by lemma y and we are left with a sum over a tree 
that we perform much in the same way we did in lemma ^. Indeed one can check that spatial decay 
appears through factors of the form 

J2 VAi^.yDA,f,XAyRXAlUA,f,XA-DA^^VAt.,, (84) 

A^,A^,A^ 

thus we can extract decay in d\{li.y,li.x) time a bound in JJ ||^a,^Va;. -C'a,^|| when we combine all 
these factors. 

Yet we need some extra features to deal with the product of Oj's, each of them being bounded in 
norm by 0(1) /x"^||Z?A,^yjZ?A,^||||Z?A,AiVfc.Z?A,Ai|| for some ki. 

The factor /i~^ can be controled with a small fraction of the probabilistic factor attached to the 
cube Aj 

If a given Df^^^ViDp^^^ appears at a large power it has necessarily a large number of links attached 
to it. Because of the tree structure, the links must go further and further so that the decay of the 
links together with the gaussian measure allow to control the factorial coming from the accumulation 



of fields. This is quite standard and the reader can refer to |Riv2] for instance. 



Finally we can write Tout.in as a sum over polymers of the form 

Tout,in = XAoutiUA,f, + 5T)xA,„ (85) 

^T = ^,^r-iH^l A ^^ns E A-i^iryT(y) (86) 

[l+L dA{Ain,Aout)] %e^,(A„„A„„0 

where ci and C2 are small constant, Fy has decay in the spatial extension of Y and 
II Yl r(y) II is bounded. ■ 

Y&A*{Ai„,Aout) 

4 Anderson model with an infra-red cut-ofF in dimension d = 2 

We are interested now in the particular case 

H = -\ + \riEVr]E (87) 
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where 



Aj7^ is a ultra-violet regularized inverse Laplacian, that we will note —p^ 

rjE is an infra-red cut-off that forces |p^ — E\'^ A}?\ logAp 

V has covariance ^ which is a Cq° approximation of a (5-function 

M^' ^ is an even integer greater than 2, and Jq G IN is such that 

M^^o s; inf |p2 - £:| <; M^^^o^^) 

Supp(r?£;) 



(88) 



For each ^ j ^ ioi we construct a smooth partition of unity into cubes of side M^ which form 
a lattice JDj. It follows a decomposition of V in fields Va and we will assume for simplicity that for 
j < k and A^ £ Dk 

V^,= Y. ^A, (89) 



even if it is not totally true because of irrelevant border effects. 

4.1 The matrix model 

We make a partition of unity according to the size oi p^ — E thanks to a function fj which satisfies 

• r/ is in Co°(]R-|-) with value in [0, 1] 

• f} has its support inside [0, 2] and is equal to 1 on [0, 1] 

• the C°° norm of the derivatives of r) does not grow too fast 
Then we construct 



%(p) 
fjj{p) 



l-ry[M2(p2_£;)2] 



fi[M^^{p'^ 



Ef] - fj M20+i)(p2 _ ^)2 for j > 



(90) 



In order to shorten expressions, we assume that 



VE 



JO 






(91) 



We expect that most of the physics will come from the neighborhood of the singularity p^ = E 
of the free propagator. As an operator in momentum space, V has a kernel V{p,q) = y{p — q)- But 
since p and q have more or less the same norm, there are only two configurations which give the sum 
p-q (cf IIFMRT2II ). 




^ k=p-' 
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We can see this in another way. If we make perturbations and integrate on V we will get Feynman 
graphs with four-legged vertices where the incoming momenta have a fixed norm and must add to zero 
(or almost zero) because of (approximate) translation invariance. Then the four momenta approxi- 
mately form a rhombus which happens to be a parallelogram. It implies that they must be more or 
less opposite 2 by 2. Thus the problem looks like a vectorial model because the angular direction of 
the momentum is preserved. 

In order to have this feature more explicit, we decompose the slice S-' = Supp(f/j) into M^'"^ 
angular sectors. We introduce fis with 

• fjs is an even function in Cq°(IR) with value in [0, 1] 

• fjs has its support inside [— 1 — -^, 1 + -^] and is equal to 1 on [—1, 1] 

• -nsil + x) = 1 - r) (l + ^ - x) for |x| sC j^ 

• the C°° norm of the derivatives of rjs does not grow too fast 

Then we define 9j = vrM"-''^ and construct sectors S^, of angular width %(! + jj) centered around 
ka = e*" (identifying IR^ and(D), with a G dj Z/, ,jy2-^. 

V, = T.(€f where [f^f)' [\p\e^^) ^ fj^) Vs{^] (92) 

Afterwards, we define the operators ?7^'s by their kernel 

Vi{x,y)=ldpe^^^^'y^fii{p) (93) 

They form a positive, self-adjoint partition of identity. 

^ = EW (94) 

We will map our problem to an operator-valued matrix problem with the following lemma whose 
proof is quite obvious. 

Lemma 7 

Let 7i be an Hilbert space and suppose that we have a set of indices I and a partition of unity 

I = T.vf (95) 

where I is the identity in C{7i) and the rji 's are self-adjoint positive operators. 
For all i & I, we define 

m = 7]i{n) (96) 

then Tt and C{7i) are naturally isomorphic to (x)'Hj and C (x)'Hi thanks to 

i(^X Kiel I 

X G W ^^ {xi)i£x where xi = r/jX (97) 

A £ C(7i) ^ {Aij)ij!zi where Aij = ijiAijj (98) 
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In our case, we define Xj as the set of sectors in the slice j and X = UXj so that we can construct 



the operator- valued matrices Va's as 



(VA)i'^ = rf^V^4 



For a slice S', we define the enlarged slice 



t} =[}T/ 



(99) 



(100) 



m>l 



Then an angular sector 5"^ of S' has a natural extension into an angular sector S*^ of S' and we have 
the corresponding operator f/^,. 



jk^ 



4.2 Size of the V^^'s 
Let V-^ be defined by 



V- 



jk 



Im 



afi 



I T/ ^m 



SjlSkmVa^Ar]}!) 



where the 6 are Kronecker's ones. We can remark that 

{vty=v^' ^ \\vi'\\ = \K'\\ 

Then we have the following large deviation result 

Theorem 3 

There exists constants C and Ce(e) such that for all A, j ^ k, a ^ 1 and A G E)k 

where C^ behaves like 1/e. 

Proof 

We use the bound 



I * A II *== -"^^ 



(vi^) (vi^ 



mo 



where 



, Im 

^13 



\ ml 
^/3q 



(At)"" 

TrA = Y.^,Al, = Y,JK^{x,x)dx 



Thus for any ttiq 

^ ' "' II * A 11?='*' 



CM-i/^ 



dH^{V) 



€ 



< 



1 



(aCMi/2)2™o 

1 

(aCM-i/2)2'no 



ivil^^^d^^jy) 



Tr 



Y'^^ M 



dH^{y) 



Let us note 



T = Tr 

-'-mo — 



vi' 



K)* 



tr 



-1 mo 



a /3 



We have the following lemma 



(101) 
(102) 



(103) 



(104) 

(105) 
(106) 

(107) 

(108) 
(109) 

(110) 
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Lemma 8 

There exists a constant C such that for all mo we have the following bound 



<Imo>^ (^2mojv^-imo 



l + M-""o(2-§)mo! 



(Ill) 



This lemma is the core of the demonstration but its proof is quite long so that we postpone it until 
the end of this part. It leads to 



P (llvii ^ aCM-^/^) ^ a-2'"o [l + M-^^^^t-i^mo! 



rfc_2 



We take mo = a M^ 3 and use the rough bound 



! < n^'^^^^e^^'^^^^ 



n\ ^ n 



:ii2) 



(113) 



to get the desired estimate. ■ 

In fact, in the proof of lemma ^ it is easy to see that rjj and % can be replaced by ry^ and % with 
the same result. Furthermore, thanks to the locality of V and to the decay of the 77's , the sum of 
several V^ 's is more or less an orthogonal sum. More precisely, for any cube Aq we define Dm(Ao) as 
the set of cubes of JDm which are contained in Aq. Then given two sets ili and 0,2 and their smoothed 
characteristic functions xui ^^^ Xn2 ^^ have 

Lemma 9 

For any n and C there is a constant Cn such that for any j ^ k and Aq G SD^. 



|Xni??iVAo%Xn2 



< 



[1 + M-"id(17i, Ao)"] [1 + M-"'=(i(f^2, Ao)"] 



max 



|%V^Ao%l|,sup sup M-^(^-"^+'=-")||r?^yAr?r, 



™</ AeD™A„(Ao) 

n<fc 



(114) 



where m f\n = min(m, n) . 

Proof 

We introduce Xaq ^ ^0° function equal to 1 on the support of Vaq then we write 

XnifijVAomXn2 = Xni%XAo%^Ao%XAo%Xf72 + Yl XniVjXAoVmVAoVnXAo^kXQi 

m<j 
n<,k 



(115) 



Afterwards we introduce the sectors and the matrix formulation and we notice that when we want 
to compute for instance the norm of the function ??"XAo^aXri2; momentum conservation tells us that 
we can convolve Xa^ ^y ^ function which is restricted in momentum space to the neighborhood of 
SJ^ — S^. In this way it is quite easy to see that we can extract at the same time spatial decay and 
momentum conservation decay. ■ 

4.3 Proof of theorem Q 

Let Aq G ^jo, we call Xc^^a and Yc^^a the events 

^c.,a = 3i ^k,3AG Dfc(Ao) s. t. ||%Fa%|| > aC.^M-^/'^M^ 
Ycy,a = \\\DA,,^VEVAor|EDA^^\\ ^ aCyjoM^''/^ 






(116) 
(117) 



19 



We will note Z the contrary event of Z. 
Theorem H tells us that 






n^c^a) ^ Y.Y. E c'e-^^*^^^^^"^^ (118) 

k j^fcAeDfcCAo) 
^ ^0(l)M20o-'=)e-t"'^^^A/^^ (119) 



1 9 ^0 

^ C7ie-2" *^^ (120) 



One can see that thanks to lemma P, Xc^a implies ^(i)c,a- Thus if we call Cq = 0(\)C 

]P{Yc„a) ^ P{Xc,a) < Cie-^-^'*^^ (121) 

Furthermore, if we work with respect to Xc^a which is stronger than Ico.a everything goes as if 
one had 

||XA,I)A,,l/A.I)A,,XA3ll ^ [1 + L-ldA(Ai, A2)]"^ [1 + L-H^{^,, Ag)]"^ ^^^^^ 

Thus we will be able to apply theorem || with an effective coupling constant 

Aeff = AjoM^o/2 (123) 

and a length scale L = M^". 

If we want to make perturbations it is clever to perturb around the expected Green's function 
without cut-off, i.e. we write 

^ ^ (124) 



p"^ - E - ifi + XV p^ - E - ino + \V + iSfi 
where //q is the expected contribution of the tadpole given by the self-consistent condition 

fio = A^Im / ^ i —dp (125) 

J P^ - E - IflQ 

Afterwards, when we compute the perturbative expansion, the tadpole with cut-off will eat up a 
fraction X^M^ ~ 0{\ log A|^^) of the counter-term so that 

^ ^ p^-E-ir]EO{X^\ log X\-^)r]E ^^^^^ 

■ 
In fact since the tadpole has a real part, it implies that we should also renormalize the energy by 
a shift 

dE = (a^ log[UV cut-off scale]) (127) 

4.4 Proof of lemma ^ 

We will note Jq = (ry^)^, Kf^ = {rfgY and X as either J oi K. 

We can perform the integration on Va so that (Imo) appears as a sum of Feynman graphs. 
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(^™o)= E /^V-<^'^' Xr^\ =E'^(^) (128) 




^■■■"'"o / IR... I \ G 



where a solid line stands for a Jq,. , a dashed line stands for a X/j. and a wavy line represents the 
insertion of a Va- In the following, we will prove the theorem in infinite volume with V having a 
covariance 6 in order to have shorter expressions. The proof can then easily be extended to short 
range covariances and finite volume except for the first few slices where one must pay attention to the 
ultra-violet cut-off but this is irrelevant because it will cost only a factor 0(1). 

The integration on Va consists in contracting the wavy lines together, then both ends are identified 
and bear an extra XA which restricts their position. 

The X's will stand as propagators and the contraction of the Va's will give birth to 4-legged 
vertices. 

4.4.1 Momentum conservation at vertices 

First, we notice that if we note a the opposite sector of a 

X^{x,y) = Xs,{y,x) (129) 

Then we put an orientation on each propagator, so that if a X^ goes from a vertex at z to a 
vertex at z' it gives a Xa{z,z') = Xa{z',z), i.e. it is equivalent to have an incoming Xq. at z and 
an incoming X^ at z' . Now, for a given vertex with incoming propagators X^j^, X^j, X^,^, X^^, the 
spatial integration over its position gives a term of the form 



^cxi...cxAx1,X2,X^,Xa) = Xa^{xi, z)Xa2ix2, z)Xa:i{x3, z)Xaiix4^, z)xAiz) dz (130) 

In momentum space, it becomes 

rai...a4(pi,...,P4) = Xa^{pi) . . . Xa^{p4) XA{k) d{pi + . . . + P4 - k) dk (131) 

where we use the same notation for a function and its Fourier transform. 

In x-space, XA is a Cq° function with support inside a box of side 0{1)M^", it means that in 
momentum space, it is a C°° function with fast decay over a scale M^-'". Thus for all n there exists 
C^{n) such that 



Cx(n)M2jo 
(1 + jVf2io|i.|2^"+i 



\xAik)\ ^ ,^ ::.:.,,, o.n+i (132) 

We make a decomposition of XA 



io/2 
XA{k) = E X.(^) (133) 

where xo has its support inside the ball of radius 2MM^-'", Xjo/2(= Xoo) has its support outside 
the bah of radius M^o/^M"^" and Xs forces |A;| to be in the interval [M^'M"^"; 2MM^M~^o]. In this 
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way, we can decompose each vertex v into a sum of vertices Vg, where a vertex Vg forces momentum 
conservation up to 0(1)M*M^-"' and has a factor coming from 

\Xs{x)\ < C;^(n)M-"" X M-*" (134) 

We split the factor in order to have a small factor per vertex and yet retain some decay to perform 
the sum on s. 

4.4.2 Tadpole elimination 

A graph will present tadpoles when two neighboring 1^'s contract together thus yielding a X{z,z). 
Suppose that we have a j-tadpole then at the corresponding vertex we will have something of the form 

dzxu{z) Kfs{x, z) Ja{z, z) Kf3>{z,y) (135) 

Between the two K^s, momentum will be preserved up to 2M M^M~^^ which in most case is much 
smaller than M~-^0' ^ so that /?' is very close to (3. Then we would like to forget about Jq, by summing 
over a and see the whole thing as a kind of new Kp. Now if per chance the new Kp makes a tadpole 
we will erase it, and so on recursively. 

First, we define the propagators as propagators (or links) of order 

(136) 
(137) 



Ot(0,0), 



Kf3^{zi,Z2)... y^^l^^{zp)xUp{zp)Ki3,{zp,y) (138) 

where we don't write the momentum conservation indices for shortness. We have a similar definition 
for ^^'Jf^a' (obtained by erasing q K-tadpoles of order 0). 

We will note X"'* to indicate that momentum is preserved up to tM~^^' between the leftmost and 
the rightmost X^s or X^>'°° if momentum conservation is worse than M"^'^''^. 

Now, we can iterate the process in an obvious way. Yet, we must add an important restriction: we 
will erase a X"'°° tadpole only if it is attached to a foo vertex. 

Lemma 10 

There exist constants Ci and C2 (independent of j and k) such that for any tadpole obtained by erasing 
a total of p J -tadpoles and q K -tadpoles we have the following bond 





Ot(0,0) _ 


Ocia' J OL 






^ 13/3' K 13 


Then we define links of order 1 






/3i.../3p- 


ai...ap J 


...dzpR 



X(;f)(z, z)\ ^ CiCf+'^M-P^-'?'=M-3-/2^(X), ^=\i'lfx = K ^^^^) 



where J^{X) is a small factor coming from the various Xs that appear in the expression of X. Thus, 

J-' gets smaller as momentum conservation gets worse. 

Proof 

First we will prove this result when momentum is well preserved, i.e. up to M^-"'" at worst, then we 

will see what has to be adapted when there is a bad momentum conservation. 
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The proof is by induction on the order of the tadpole. We define Ci, C2 and C3 such that 

\X^{x,y)\ sc C7iM-3^/2 ^140) 

sup/dy |X^(x,y)| ^ C3 (141) 

C2 = 9C1C3 (142) 
It is easy to see that for level tadpoles 

^ CiM-3^/2j^(X) (143) 



Jx^^i \Z, ZJ 



Now, consider ™J^^/ a J-tadpole of order m and weight {p, q) obtained by erasing n i^-tadpoles 
of order m — 1 and weights {pi,qi), ■ ■ ■ , (pn, Qn)- We have 

p = pi + ...+Pn q = qi + ... + q„ + n (144) 

The expression of "J will be of the form 

"'J!^\^,^) = T. T. jdz^...dz^Uz,z^)^-^-^Kf^^^\z^,z^) 

Ql...Q„-l fii---fin 

XuA^l) Jai izi,Z2) ■ . . K^p^^^"\zn, Z^) X«„ (^n)^' (^n, Z) (145) 

Since we supposed that we have momentum conservation up to M"-'"'^, the Qj's will be either 
aj_i or one of its neighbors and /3^ will be either /3j or one of its neighbors. Thus the sum on sector 
attribution will give a factor 32™- Im"*^/^ ^ 9"M"^/^ 

We have n + 1 J's but only n spatial integrations because we have a tadpole. This gives a factor 
C'^CiM^'^^i'^ (we forget about the momentum conservation factor for the moment). 

Finally the ""-^i^'s bring their factor so that 

^ Ci(9CiC3)SP'+S'?'+"M--'Tp>-^'(E'?'+«)M-3i/2^(X) (146) 

which is precisely what we want. Then we can do the same for the ™i^'s. 

Now we must consider the cases with bad momentum conservation. First, let us suppose that 
momentum conservation is bad overall for '"J but was good for the ™~^i^'s, then the previous argument 
will work except if there are some v^ vertices. In this case we will have to pay a factor M^'"^ to find 
the following Qj instead of a factor 3. But from the corresponding Xoo we have a small factor 

1 + MioiV'/2 (^^^) 

from which we can take a fraction to pay the M^'"^ and retain a small factor for J^{X). 

Finally, if a "^^^K has a bad momentum conservation it is necessarily attached to a Voo vertex 
(otherwise we would not erase it). In this case we must pay a factor M^''^M ''^ (to find f3'^ and ctj) 
but again we can take a fraction of the factor of Xoo to do so. 

When tadpole elimination has been completed, we have erased tj J-tadpoles and t^ i^-tadpoles 
and we are left with tjiq = mo — tj — t^ vertices linked together by ttiq J's and mg K's (a tadpole 
which has not been erased being seen as a propagator). 

For a X^^r{x,y), it is quite easy to see that to integrate on y with fixed x amounts more or less 
to the same problem for 0{iy~^'^Xa and that to find a' knowing a costs a factor 0(1)^"'"'^. 

23 



4.4.3 Sector conservation at the vertex 

Lemma 11 

Let (S^^^, . . . , Sl^^) be a quadruplet of sectors of the enlarged slice S^ and ^ r < 0(1)M''^ such that 
there are pi G S*^^ , . . . , p4 G 5^^ verifying 

\pi + .■■+P4\ ^ rM'^ with (148) 

Then we can find {q, a', /3, (3'} = {ai, . . . , 04} satisfying 

{ \a'-a\ ^ ia^ + b)M'^/^ 

\ 1/3' -^1 ^ (aVF + 6)M-'/2 ^'^^^ 

where a and b are some constants independent of I and r. 

Proof 

If we can prove the result for I ^ 0(1) then we wih be able to enlarge the result to any / provided 
maybe we take some slightly bigger a and b. Therefore we assume that this is the case in the following. 
We define {a,a',P,f3') by 

• {a, a', (3, (3'} = {ai,...,a4} 

• a = ai 

• \a — (3\ = min la — aJ 

JG{2,3,4} 

• \a — a'\ ^ |a — /3'| 

Then, if |a — /3| < \a' — /?'| we exchange (a,/3) and (a',/3'). 

A sector S is included in a tube, of center ky = e'"'^ and whose direction is orthogonal to the 
direction 7, of size 



length : L = tiM-^I'^(\^ ^) 



\ width : 2M-^ 
We define 



(150) 



k^p = kc + kf3 = 2cos[ - — - ) e'^ = 2 cos X e^^ = re'^ (151) 



k^if^i = -ka'i3' ^2 cos X e (152) 

If we can prove that 



\x' - x\ sC (a'^+5')M"'/2 
\e'-e\ ^ {a'^ + b')M~^/^ 



(153) 



then we will be able to conclude, with a = 2a' and b = 2b'. 
It is easy to check that by construction, we have 

• < x' s; X 

• |a - /3| ^ ^ ^ cos X ^ i 
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We have a trivial bound 



^0/3 — f^a'/S 



^r\' R'\ ^ ^ tan 



s$ rM'^ + 2L + 4M'^ = R 



Therefore 



W -e\ <2tan 



^R^O{l)M 



-1/2 



(154) 



(155) 



We can see that is very well conserved. 

If sin X < (ai^ + 6i)M-'/2 then |x - x'| s$ x ^ (02^/^ + 62)M"'/2^ 

Otherwise, let us remark that S^ + 5*1 is at a distance at most 2M^' from a rhombus Rap of center 
kaj3 and of diagonals 



2L sin X in the direction u^ = ^^-^ 

2L cos X in the direction ug = ^^^ + ^ 



(156) 



Then, Rap — Ra'B' i^ ^^ ^ distance at most AM from a rectangle 7Z of center /c^^ — k^igi and of 



sides 



Lr = 2L (sinx + sin x' cos \9' — 9\ + cos x' sin \9' — 9\) in u^ 
Lg = 2L (cosx + cos x' cos \9' — 9\ + sinx' sin \9' — 9\) in u^i 



(157) 



Since \9' — 9\ < 0{1)M '/^, we have | cos{9' — 9) — l\ < 0(1)M '. We define a z axis in the direction 
Ur- {ka'pi — kajs) has a z coordinate 2(cosx' — cosx) + 0{1)M^K This leads to the condition 

2| cosx' - cos x| < rM~' + 2L(sinx + sinx') + fesM"' 
< (r + 63)M"' + 4Lsinx 

Let us note that | cos(x — u) — cosx| is an increasing function of u and that we have 

cos(x — u) — cosx = sinxu cosxu^ + u^eiu) with \e{u)\ < - 

2 6 

We take u = {^ + bi)M-^''^ + 2L = {^ + b'^)M-^/^ < 0(M-'/^). 

I / \ I ■ /I ""A 2 

cosx — u— cosx ^ smxu— — I — \u 

^ 2L sin X + sin x( Vf + bi)M-^/'^ - {\/r + b'^fW^ 
Since we are in the case sinx ^ (01-^? + bi)M ''^, for ai and 61 large enough we will have 

2| cos(x — u) — cos x| ^ (r + 63)M~ + 4L sin x 
Therefore we must have |x' — x| < u < (v^ + b'^M^^i"^ which allows us to conclude. ■ 

4.4.4 Size of a graph 

The previous section shows that, at each vertex, momenta come approximately by pairs of opposite 
sectors. Thus, for all the vertices which haven't been erased by the tadpole elimination process, we can 
choose by a factor 3 how to pair the sectors. Then we split the vertices in two half- vertices according 
to this pairing. We represent graphically this as 



(158) 
(159) 



(160) 

(161) 
(162) 

(163) 
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a p 

P P 

This gives 3™'o (split) graphs that we will consider as our basic graphs in the following. 

A graph is decomposed into a number of momentum cycles connected together by wavy lines. We 
will follow those cycles to fix momentum sectors. Finding the enlarged sectors (of level j) will cost a 
factor M^'"^ per cycle times a constant per vertex. Then we will pay an extra M^ ^^''"^ for each K 
propagator to find its sector. 

We define c the total number of momentum cycles that we decompose into t tadpoles, b bubbles 
(with 2 vertices) and / large cycles (with 3 or more vertices). We have 

t + h + l = c (164) 

i + 26 + 3/ s; 2mo (165) 

and the sector attribution costs 

Ai = C"'0M^j/2^jm',{k-j)/2 ^^gg) 

Notice that the constant has an exponent niQ because of the tadpole elimination process. 

The spatial integration of the vertices will be made with the short J links whenever possible. We 
can decompose each graph into J-cycles linked together by K links (because there are 2 incoming J's 
at each vertex), this allows to integrate all the vertices but one per cycle with a J link. The total cost 
is (noticing that the last vertex is integrated in the whole cube A) 

A2 = C™oM^('"o"'^')-^'/^M^('='~^)*=/^M^'' (167) 

where c' is the total number of short cycles that we decompose into t' short tadpoles, b' short bubbles 
and /' large short cycles. 

t' + b' + l' = d (168) 

t' + 26' + 3r sc mf) (169) 

The scaling of the tadpoles and the propagators give a factor 

Tadpoles that have been obtained by erasing a few vertices (say 0(M-"''^) for instance) will have 
an extra small factor because they strongly violate momentum conservation, we can take it to be a 
power of M~-^0'^. Tadpoles with higher weights will not have this good factor but we will see that 
they bring a better combinatoric. The t momentum tadpoles will consist in t\ low weight ones and t^ 
others while the t' short tadpoles split into t'^ low weight ones and t'^ others. We can manage to have 
a factor 

If we have a short bubble we will have four incoming long propagators whose momenta must add 



up to zero up to xM~^'^. If we apply lemma 11, we can see that knowing 3 of these momenta it cost 
only a factor 0{l)^/x to find the fourth momentum sparing us a factor M^ ^^''"^ obtained by naively 
fixing first the enlarged sector at slice j. If the bubble has a weight p, i.e. the two short propagators 
have been obtained after erasing p vertices, and a momentum conservation worse than 0{l)pM~^° 
then the small factor of bad momentum conservation will pay for the 0{\)y/p. We will have b'-^ such 
good bubbles, each of them bringing a factor M ~' ■'^' ^. In addition, we will have 63 bad bubble of 

26 



weight greater than M^^ ^' for which we earn nothing and 63 bad bubbles of weight pi bringing a 
factor C y/plM^^^'^^'i'^ . This gives a factor 

A^ = M-(''i+''3)('=-i)/2c''3 [] ^ (172) 

i 

Finahy we have the following bound for the contribution of a graph 
\A{g)\ ^ Ai...A^ 

]^-3m'giJ+k)/2j^-2jo(ti+t',)]^j-{b[+b'2)(k-j)/2 TT ^ (^-^73^ 

i 
^ C"^0 ]^k/2 jY^cj /2 j^3c' (k-j) /2 j^j-in[^{k+j /2) j^-jtj-ktk) j^~2jo(ti+t[) 



If we use equations (165) and (|169|) we obtain 



Mt2J/3Mt'^{k-j)^jb'^{k-j)/2^-t^{2jo~j/3)^^t[{2jo+j-k) -Q ^ (175) 

i 

we will take mo ^ M^'^ so that M^'"^ < C™". Furthermore, t and t' are at most equal to M~-^0'^m,o 
thus M^* ^ 6""°. It allows us to rewrite the bound 

i 

4.5 Graph counting 

Lemma 12 

-Lei T(p) 6e the number of ways to contract 2p adjacent V 's so as to make only generalized tadpoles. 
We have 

Proof 

It is easy to see that a good contraction scheme, i.e one that gives only generalized tadpoles, corre- 
sponds to have no crossing contractions. It means that if we label the fields Vq . . . V2p-i according to 
their order and if Vi and Vj contract respectively to V^ and Vi then 

i<j^k<j oi k>l (178) 

We have r(l) = 1. For p > 1, we contract first Vq to some Vi. Vi, . . . , Vi^i will necessarily contract 

among themselves making only generalized tadpoles and so will do V^+i , . . . , V2p-i . Thus i is necessarily 

odd and we have 

p-i 

Tip) = J2 T{k)T{p -l-k) (179) 

A:=0 

where by convention r(0) = 1. 
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We introduce the generating function 

oo 

tiz) = Y^T{p)zP (180) 

p=0 

The recursion formula ( |179| ) can be translated into an equation for t which is 

t{z) = zt^{z) + 1 (181) 

whose resolution yields 

i + vr^4i i-yr^rii 

t{z) = or (182) 

Since the second solution is analytic around z = 0, we can take it as t{z) and the coefficients of its 
power expansion will give us T{p). An easy computation leads to the desired formula. ■ 

Lemma 13 

The number Nm{B) of graphs with B possible momentum bubbles obtained in the contraction of a 
cycle of 2M V 's has the following bound 

Mm{B) ^ C^'{M - B)\ (183) 

Proof 

First lets us remark that bubbles come in chains (possibly with a tadpole at one end) of two possible 

types 

• typ- 1^ xnxix 



tyP« 2: X±X±1J^ 



where a solid line stands here either for a J or a iC. 
We have two special cases 






* ^ X A A X y which can be seen as a type 1 chain 

• ~\/j~ which can generate only one momentum bubble so that we can see it as a type 2 chain 
of length 1. 

Having chosen the V's, there are only two contraction schemes that yield a type 1 chain and a 
unique contraction scheme for a type 2 chain. If we fix explicitly the subgraphs corresponding to B 
bubbles and contract the remaining l/'s in any way we will get all the desired graphs plus some extra 
ones so that we can bound Mm^B). 

We construct ri type 1 chains of lengths /3i , . . . , /3ri and r2 type 2 chains of lengths 71 , . . . , 7r.2 . 
We set 

Bi = J2/3i (184) 

i 

B2 = J2^i (185) 

i 

B = B1 + B2 (186) 
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To count the contraction schemes, first we cut the cycle of 2M V^s into a sequence of F's (there 
are 2M ways to do so). Then it is easy to check that in order to build a type 1 chain we must choose 
two sets Bi and 8i of /?j + 1 adjacent V^s while for a type 2 we need a set Vi of 2^i + 1 adjacent V^s. 
We distribute those 2ri + r2 objects in (2M — 2Bi — 2ri — 2B2 — r2) + 2ri + r2 boxes in an ordered way, 
and for the i type 1 chain the respective order of Bi and Bi will fix the contraction scheme. Then, 
there remain 2M — 2B — 2ri V^s to contract so that we have the following number of configurations 

U„[B),,2M Y. T. Y. ;^;^(2ra^^^r^P*'-2B-2n-l)!! (187) 

B,+B2=Bri<Si ,3i + ...0,j=SiAS»i '^i-^2- l-^JVJ ^o zri r2). 
We can compute this 

Mm{B) < 2^^^^ 1^^^ ^^^ l^r, - 1 j l^r, - 1 j [2(M - B) - 2r, - r,]!(2r,)!r.! 

'■2<S2 

(2ri)! \2{M-B-rx)\\ 



n! 2^-^-<-^{M - B - rx)\ 



(188) 



2^^ E E f^'_ ^)f^'_^)32(^^-^)22'^Vi!2*^-^-^i(M-i?-ri)!(189) 

'■2<S2 

< 2M J2 2^i-i2^2-i9A/-B2^^(M - S)! (190) 

Bi+B2=B 

^ 2M(S + 1)18*^(M-B)! (191) 



4.6 Bounds 

Now we can achieve the proof of lemma |8| in bounding 

g 

In order to compute this sum, we fix first tj, t^ and 6, where h is the number of possible momentum 
bubbles and therefore is greater than h. 

Then we define the set Q{tj,tk, b,n,qi, . . . , qn) has the set of graphs with the corresponding tj, t^ 
and b and for which the erased tadpoles form n sets of 2qi adjacent V^s. We can write 

To bound |^(^)|/ n(9i + 1) 'we notice that when a graph has a bad bubble of weight pi it means that 
we have erased two set of generalized tadpoles qi- and 52^ on the two propagators of the bubble with 
Qii ~\~Q2i = Pi- Thus we have a corresponding factor [qi. + l)~^(g2i + 1)^^ which control the bad factor 
^Jpl of the bad bubble so that 

h+tk , 
<^mQ>^ Y. H —\ Yl E \{{qi + l)C'^''M-<^^+^^/'^M~^^^~^*'^M-'^<~^'^^'^ (194) 



n 

1' 



t,,t^j, n=l ■"• 9i,.^9n ggo(..,) 
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The number of graphs in i}(. . .) has the following bound 
This leads to 



tj .t-L. 



E;^ E 1 (196) 



n! 



n = l qi + ...q„=tj+ti^ 



tj +tk — I 

n — 1 



^ ^ C"«M-"o(^'+i)/2M-i*.-^*^M-K-%/6(^^_6)! ^ ''^■^'''^ (197) 

^ J2 C"oM-™o^M-J*^-'=*^M-K-b)§(m^-6)! (198) 



^7 5^£- h 



Summing on tj and tk is equivalent to sum over ttiq and tfc with tj = mo — tfc — JTig. The sum over 
b is roughly evaluated by taking the supremum over b, the result depends whether ttiq is greater than 
MJ/^ or not. 

+ E E'^"'°^~'""^'''^-'^^^^"-'*'"''*'™ax[l,M-™oi/6(^[j_5)!l (199) 
sc C'"oM-'"o^' ^ j^-mi(fc-i)/2^-tfc(fe-i) r-^ ^ 1^^, ^^,./6^M-'"oi/6(^^ _ ^),1 (200) 



m'r,,t^ 



Finally, the sum over i^ is easy and we bound the sum over itiq by finding the supremum. One 
can check that it gives the announced result. ■ 

5 Conclusion 

Understanding the effect of perturbations on the free spectrum of Hamiltonian operators is an out- 
standing challenge. We think that for the two-dimensional case, this paper can help to control the 
model up to quite close to the scale of the expected "mass" {i.e. imaginary part) so that in studying 
the full model, one can focus on the thin slice of momentum p^ — E ^ }?. Then we think that a key 
of the problem lies in the fact that in this case the potential is very close (in momentum space) to 
a large hermitian random matrix with almost independent entries and therefore we can connect our 
problem to the much better understood domain of random matrices or equivalently use the vector 
model picture of the problem. 

One can note that in dimension d = 3, the interpretation in term of random matrices can also be 
helpful (cf. [MPR| ) , but then one has to deal with constrained matrices whose entries are no longer 
independent. 
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